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We develop dynamical non-Markovian description of quan- 
tum computing in weak coupling limit, in lowest order ap- 
proximation. We show that the long-range memory of quan- 
tum reservoir (such as the one exhibited by electromag- 
netic vacuum) produces strong interrelation between struc- 
ture of noise and quantum algorithm, implying nonlocal at- 
tacks of noise. This shows that the implicit assumption of 
quantum error correction theory - independence of noise and 
self-dynamics ~ fails in long time regimes. We also use our 
approach to present pure decoherence and decoherence accom- 
panied by dissipation in terms of spectral density of reservoir. 
The so-called dynamical decoupling method is discussed in this 
context. Finally, we propose minimal decoherence model, in 
which the only source of decoherence is vacuum. We optimize 
fidelity of quantum information processing under the trade-off 
between speed of gate and strength of decoherence. 

Pacs Numbers: 03.65.-w 



I. INTRODUCTION 

In spite of many remarkable results on decoherence in 
open quantum systems Q the problem still remains open 
involving some conceptual and interpretational difficul- 
ties in the description of dynamical quantum effect. Re- 
cently the problem become crucial in connection with the 
idea of quantum computers The latter stimulated 

a huge theoretical and experimental effort to control evo- 
lution of quantum systems. This progress allows to hope 
that theoretical obstacles to build quantum computer 
(QC) can be overcome. To deal with unavoidable de- 
coherence, the quantum error correcting codes have been 
designed (|J2|] resulting in the theory of fault-tolerant 
(FT) quantum computation In this theory it is 

tacitly assumed that decoherence acts independently of 
the structure of the controlled self-evolution of QC (the 
one including quantum algorithm and the scheme of er- 
ror correction). The self-evolution only propagates the 
errors. In fact, decoherence depends on the kind of dy- 
namics of the system of interest, and this should be taken 
into account. 

The main purpose of this paper is to provide more 
complete general analysis of interrelations between deco- 
herence and controlled evolution, basing on well estab- 
lished theory of open systems . We develop the dy- 
namical, Hamiltonian description of decohering quantum 
computer, deriving and analysing non-Markovian master 
equation. We show that the long-range quantum mem- 



ory in conjunction with self-dynamics of quantum com- 
puter implies highly nonlocal structure of noise. This 
has remarkable implications on quantum error correction 
concept. The key idea of the latter is that the noise is 
local, so that the quantum information can be hidden in 
multi-qubit entanglement. Once we do not couple the 
qubits through quantum gates, the noise is local indeed. 
However, to protect the information against the noise, we 
need to produce multi-qubit entanglement. Due to mem- 
ory, the reservoir will then "see" the evolution as highly 
nonlocal, and after some time it will act nonlocally it- 
self. We argue that the long-range quantum memory is 
unavoidable for fundamental reasons: it is exhibited by 
interaction with vacuum, which can be never removed. It 
turns out that the memory is relevant even for quantum 
optics regime, despite the fact that the deviation from 
exponential decay is unobservable. 

In Refs. §1 it was noted that FT method may not 
work, if time or space correlations do not decay exponen- 
tially. What we show here, is that correlations in time 
cause the noise to follow the evolution of the system that 
is to protect against the noise. It follows that FT scheme 
should be revisited to take into account the memory ef- 
fect. There is a hope that it is possible, as the collective 
noise attacks due to memory are of special type: they can 
be viewed as a result of propagation backward in time of 
single-qubit errors. 

One should mention, that the long quantum memory 
is exhibited also by phonon environment which charac- 
terized stronger coupling than electromagnetic vacuum. 
This means that within the program of solid state quan- 
tum computation the memory effects can become practi- 
cally relevant. 

Other programs to avoid decoherence use dynamical 
decoupling [ pd] |, and decoherence free subspaces Q. We 
discuss the first approach in the context of real physical 
environment and conclude that it is hard to find condi- 
tions for which the method could be useful. The decoher- 
ent free subspaces approach bases on opposite paradigm 
to the one of quantum error correction - it exploits sym- 
metries of collective interaction with environment. Even 
though our approach is fully general, in this paper we 
will discuss the independent interaction case. 

We also introduce a minimal decoherence model of evo- 
lution, within the framework of non-Markovian approach. 
The motivation behind the model is to try to take into ac- 
count only fundamental obstacles for building quantum 
computer. That is we remove any decoherence mecha- 
nism that could be, in principle, avoided at some tech- 
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nology. In the model the only environment is the vac- 
uum, and there is no internal, natural self-evolution (the 
only evolution of the system is constituted by quantum 
gates). It follows that there is some room to optimize 
the process of quantum computation within this model. 
Subsequently, we obtain time-energy trade-off due to that 
fact that to perform a given gate one needs a constant 
amount of "action" that can be split in different ways 
into time and energy. Higher energy implies quick gate, 
but simultaneously, enhances decoherence. In result it is 
possible to keep high fidelity at the expense of worse scal- 
ing. For sequence of n gates, the physical time needed to 
keep high fidelity scales as n"^/^. 

This paper is organized as follows. In Sect. |l|we make 
some preliminary remarks. Then we develop Hamilto- 
nian approach in Sect. |l|. We obtain the transition map 
(called here error map). Then we pass to the Markovian 
limit for time- independent Hamiltonian and indicate that 
in non-Markovian stage the reservoir recognizes structure 
of levels of the systems. We discuss possible types of 
spectral densities causing decoherence and the problem of 
high frequency cut-offs. Subsequently we discuss the dy- 
namical decoupling method. In Sect. [V we present two 
examples - QC in memoryless reservoir, and QC driven 
by kicked dynamics. We discuss error path interference 
effect. In Sect. |IVq memory causing by interaction with 
vacuum is derived. Then we pass to the main result of the 
paper (Sect 0) der iving and discussing the master equa- 
tion for decohering QC. We relate the results to quantum 
error correction concept. In Sect. Vl we postulate 771m- 
imal decoherence model, and provide formula for fidelity 
within the model. An example of fidelity for single qubit 
rotation is presented and the time-energy trade-off is dis- 
cussed. 



II. PRELIMINARIES 

A. Fidelity in interaction picture 

We consider the general case of QC interacting with 
environment E. The initial state of QC is given by ipto- 
The dynamics of QC without decoherence, is given by 
unitary evolution of quantum algorithm 



(1) 



We will call algorithm the total controlled evolution, i.e. 
the computation algorithm itself plus the error correction 
procedures. Due to interaction with environment, the 
initial state will evolve into mixture 



(2) 



with Qtg — ipto ■ The fidelity of decohered computation is 
given by 



Ft = (HQtli't). (3) 
It is convenient to use interaction picture. Putting 

A{t,to) = U{to,t)A{t,to) (4) 

we obtain 

Ft = {i:tom,tomto){Mmo)- (5) 

Here we used the notation 

dig) = UgUl (6) 

where U is unitary operation. 

Note that in the interaction picture the quantum al- 
gorithm does not compute, but only "interacts" with er- 
rors, causing them to spread more and more. The map 
A{t,to) describes the net effect of decoherence. As we 
will see, in FT model the decoherence is "decoupled" 
from algorithm, so that the total decoherence amounts to 
the faults caused by environment, which are then prop- 
agated by the algorithm. In Hamiltonian approach, the 
interaction with environment will be "entangled" with 
algorithm, due to quantum memory, and no such simple 
description will apply. 

B. General evolution 

Through the paper we will deal with equation 
d 



dt 



at — Ct<Jt + Lcjt 



(7) 



where Ct is some time dependent operator (it will de- 
scribe self-evolution), while L is time independent oper- 
ator (interaction), at is vector. Let the free evolution of 
a be given by the operator T[t, to) 

T{t,to)at„=at. (8) 

The operator satisfies 



— r(i,to) = £tr(i,to). 

The total evolution denoted by A{t,to) satisfies 



d ~ 
di 



A{t,to) = AA(t,to) + LA(t,to). 



(9) 



(10) 



Then the evolution in interaction picture given by 
A(i, io) = ^{tQ,t)A{t,to) we have the following formal 
expansion 

A(t,io) =1+ ^ / di™... / dti X 

)L...LT{tiM) (11) 

Full evolution is given by 



A(t,to) = r(i,to)-f ^ / dt 

m=l 

r(^, tm)LT{tm, t„i^l)L . . . LT{ti, to) 



CO i,t „t2 

_ ,„.../ dt 1 X 



(12) 
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III. EVOLUTION OF QC IN SECOND ORDER 
APPROXIMATION 



In this section we derive reduced dynamics of QC in- 
teracting with environment R. Consequently, in eq. (|^) 
the operator £t will be sum of free Hamiltonians of the 
QC and R, the operator L - the interaction Hamiltonian, 
and at - the wave function of the total QC + R system. 



A. Non-Markovian Born approximation 

The Hamiltonian of the total system is of the form 

H = Hqc + Hii + \Hinf (13) 

Here Hqc is time-dependent Hamiltonian of the com- 
puter, Hr is time-independent self-Hamiltonian of the 
environment, A is coupling constant, which is assumed 
to be small (A « 1). In the following we will remove A 
from formulas by incorporating it into interaction Hamil- 
tonian Hint- The latter is of the form 



Hint — ^ Sa ® Ra 



(14) 



where Sa,Ra are self-adjoint operators. The resulting 
evolution of the total system is described by unitary 
transformation U{t, s) such that 



U{t,s)g{s)^ Q{t). 
We use the notation: 

H{g) = [H, g]. 
The self-evolution of QC is given by 



(15) 
(16) 

(17) 



where T is time ordering operator. Finally, the free evo- 
lution of the total system is 

Uo{t, s) = Uc{t, s) ® e^*^«(*-"). (18) 

The initial state of the total system is given by 

g{s) = gc{s)<S>iOR, (19) 

where gc is the initial state of quantum computer, while 
is stationary state of environment. Without loss of 
generality, one can assume that 



TrRaUJR = 0. 



(20) 



Besides, the state lor commutes with dynamics of envi- 
ronment 



We consider lowest order approximation of the evolution, 
obtaining from equation (pj 



C/(t, s) ~ f7o(i, s) - i j duUoit,u)HintUoiu,s) - 

s 

t u 

du J dwUo{t,u)HintUo{u,w)HintUoiw,s). (22) 

Then the evolution of the reduced density matrix reads 
as gc{t) = Ti R{g{t)) is given by 

t U 

Qc{t) = Uo{t,s)gc{s)-TTRi^j du j dwUo{t,u) x 

s s 



HintUo{u,w)HintUQ{w,s)gc{s) (g>UjR 



(23) 



Subsequently, we find that gc{t) can be written in the 
following compact form 

Qc{t) = Uc{t, s)(gc{s) ~ ^[A{t, s)gc{s) - gc{s)A{t, s)] 



+^t,s)gc{s)-i[h{t, s),gc{s)] 



(24) 



The completely positive superoperator ^{t, s) is given by 
the following formula 



I z 

^{t,s)gcis) J '^^ J d«^(i?/3i?a(w-u;)). 



S'/3(m, s)gc{s)Sa{w, s) 



(25) 



where 5„(m,s) ^ Uc^{u,s)Sa, Ra{t) = e'"'^* RaC'"'^* 
are the versions of Sa and Ra evolving in Heisenberg 
picture according to free evolution; (i?Q,i?/3(t))^^ — 
Ti{ujRRaRi3{t)) is the autocorrelation function of the en- 
vironment. The operator A{t, s) is given by 



A{t, s) = At(t,s) = $*(t,s)I 
where l>* is dual map to l>. Exphcitly 



(26) 



t t 



[tOR, Hr] = 



(21) 



A{t,s)^'Y J du J dw {RaRpiu- w))^jjSaiu,s)Si3{w,s). 

"/3 s s 

(27) 

Finally, h{t, s) describes the Hamiltonian contribution to 
the dynamics due to the interaction with environment 
(Lamb shift, collective Lamb shift, etc.). In the following 
we put h{s, i) = by applying renormalization proce- 
dure. We add to the Hamiltonian HQc{t) the appropri- 
ate counter-terms which cancel the contribution h{t, s) in 
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a given order of perturbation calculus. Therefore in the 
following Hqcit) is the full physical Hamiltonian con- 
taining all relevant terms. In this way, when passing to 
Markovian approximation for the reservoir at the ther- 
mal equilibrium we obtain the Gibbs state corresponding 
to full Hamiltonian, as it should be. 

One can pass to the frequency domain, putting 



+00 



Rap{uj)e-''^*duj, 



Yaiuj) = / Saiu,s)e~ 



(28) 



(29) 



The function Rap is called spectral density. We then can 
write $ as 

+ CXD 

Ht,s)gc{s)^Yl J <iujR^t,{Lo)Yp{uj)gc{s)Y^{Lo) (30) 

-00 

where we do not write explicitly the dependence on t 
and s of the operators Ya. Denoting 3>(t,0) = $4, the 
evolution in interaction picture can be thus written (for 
s = 0) as 



gc{t) = Uoit,o) 



(31) 



where go — f?c(0)- The above equation closely resembles 
the form of the generator of the Markovian semi-group 
1^,^. In next section we will exhibit Markovian ap- 
proximation of this formula. 

One notes that the crucial element is here the transi- 
tion map $ (we will call it error map) . Operators Sa are 
errors that can occur during interaction with reservoir 
(cf. jl^). In particular, the map contains propagation of 
errors which is nothing but evolution of Sa in interaction 
picture. The overall error is proportional to ||$||, where 
II • II is some norm (see B). 



B. Markovian limit for time independent 
Hamiltonian 



with 

Then the formula (BOl) reads as 



$(s,i)gc(s) 



(33) 



(34) 



Sf:i{ujk)gcis)Sl{ujk) 



TT ^ j du Rap{Lu)5f\L0 - Uk) 



S0{i^k)gc{s)Sl{ujk) 



e 2 



(35) 



Here we have used two models of Dirac delta function 
with width - 



(36) 



Approximately we have {S^^Y ~ 'kt6'Y\ The second, 
"non-resonant" term of eq. (|3^) will vanish in the limit 
of large r, as the overlap between two deltas will decrease 
for large r. (This is an alternative form of a well known 
rotating- wave approximation). More precisely, one re- 
quires T » where 

Auj = mm{\ej - ej>\;ej ^ e^} 

Finally, we obtain 

^{t,s)gc{s) ^ TTT ^ Raf3{uJk)Sp{uJk)gc{s)Sl{uJk) 

(37) 

Note that the map depends on t and s only through the 
factor T = t—s. Thus is exactly the transition map for 
the quantum dynamical semi-group in the weak coupling 
limit HITBl 



gc = Lgc 



(38) 



Let us put t = ^, s 



and 



(32) 



where {\j)} is the orthonormal basis in the Hilbert space 
of QC. To find the meaning of the completely positive 
map <I>, we will find the behaviour of evolution for long 
time. We have 



with 



Lgc = -i[Hc, gc] ~ ^{^*{^), Qc} + -^gc- 
It t 



(39) 



One can see | |l6[ that we can distinguish the initial, non- 
Markovian stage, when reservoir is "learning" the struc- 
ture of the Hamiltonian of the system. This requires 
time T » Once the structure is recognized, we 

have Markovian stage, during which the system is being 
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relaxed towards the Gibbs state of thermal equilibrium 
^Q-PH determined by H. 

Suppose now that we have two qubit system, where 
the components are coupled to one another: e.g. the self 
evolution is to produce some two qubit gate. Then, in 
Markovian stage the system relaxes to compound Gibbs 
state. Thus, the noise, after recognizing that the Hamil- 
tonian is compound, becomes compound itself We 
will analyse it in more detail in Sect. ^ 

C. Decoherence and dissipation in Markovian 
regime. 

The properties of the decoherence in Markovian limit 
can be read from the formula (|3^). We will now anal- 
yse the formula to get first intuition about decoherence. 
First of all, one can single out pure decoherence which is 
not connected with energy exchange with environment. 
The populations of energy levels are kept, but the phases 
subject randomization. The relevant term is 



t:tY,Rc.A^)SMqcSI{Q). 



(40) 



a, 13 



Indeed Sa{Q) commutes with He- Hence this term does 
not lead to transitions between energy levels. Thus the 
term i?Q/3(0) stands for the strength of pure dephasing. 
On the contrary the other terms of the sum correspond 
to the decoherence accompanied by energy exchange, be- 
cause the operators Sai^Jk) either do not exist (if Sa com- 
mute with He) or describe transitions between energy 
levels. Note that if the spectral density R^p vanishes for 
a; = 0, there is no pure decoherence: the decoherence is 
always accompanied by dissipation. Assume, for exam- 
ple, that a system does not have self-Hamiltonian (e.g. 
spin of free electron) . Then there is no dissipation. If in 
addition i?Q_^(0) = then there no decoherence at all. 



D. Important examples of reservoirs 

We note that an important characteristic of the in- 
teraction with environment is the shape of the spectral 
density. A general property of the heat bath at the tem- 
perature T is the following 



(41) 



In particular, at zero temperature, only nonnegative fre- 
quencies are relevant. Let us now present three important 
examples of environments. 

a. Linear coupling to bosonic field Operators Ra are 
given by 



Where a^(0c(), a{(j)a) are creation and anihillation op- 
erators for fields </)„, respectively. The following form of 
diagonal elements of Rap can be obtained 

Ra,a{oj)^ y"dfc5(fi(fc) - w)(n(fc) l)|0„(fe)P (43) 

where Vt{k) is the energy of the boson {VL{k) = ^ or 
Vt = vk), n{k) = ^f,nlk)_i - Typically Raa{(^) grows like 
I a; I'' for small u and then rapidly falls down for > ojc 
{ujc is a cut-off frequency). For electromagnetic inter- 
actions in dipol approximation we have Ra^a{(^) ~ to^ 
for small lu. Similarly for some cases of phonon interac- 
tion. Cut-off parameter Uc is model- dependent (see, e.g., 
I p^ ) and characterizes the range of validity of model but 
does not mean that the frequencies co > lOc are not in- 
fluenced by noise. For example, the electromagnetic in- 
teraction can be described by the dipol approximation 
for — — Ikl « — where rn = — ^-r is classical electron 
radius . For larger frequencies different model involv- 
ing non-bounded (scattering) states of electrons should 
be taken into account. The other cut-off for free elec- 
trons LOc — 2meC^ means only that for larger frequencies 
the pair production should be taken into account. In solid 
state LOc is a Debye frequency. For higher frequencies the 
coupling to phonons makes no sense and we have to con- 
sider local interactions with the neighbouring atoms or 
ions and electrons. 

b. Interaction with a dilute gas. We use the model of 
free bosonic or fermionic gas in the low density approxi- 
mation 

Ra = a\4>a)a{(l3a) - (a^((/'a)a(0a)) (44) 
The spectral density is then of the form 

Raaiuj)^ J j dkdln{k)\(j)a{k)\^\Ml)fSi^ik) - ^(l) - 

(45) 

Typically Raa(0) > and, again Raai^j) falls down for 
> ujc - & model dependent cut-off frequency. 

c. Fluctuations of molecular field. Influence of a lo- 
cal fluctuating field can be described by a classical noise 
Fa(t). For example for a coloured noise model we have 



(FMt)) = De-- 
Then the spectral density is Lorentzian 

D 

Raa{^) — 



(46) 



(47) 



Ra = a{(j)a) + «■'■((/)„) 



Such a density is relevant for relaxation mechanism in 
pulsed NMR experiments. 



(42) 



5 



E. Discussion of dynamical decoupling method 

From eq. we see that one can gamble with deco- 
herence by changing spectrum of the Hamiltonian trying 
to choose frequencies LUk for which the spectral density 
is small. In Ref. pl the method of dynamical decou- 
pling was exhibited. It bases on adding periodic, rapidly 
alternating term to the self-Hamiltonian of QC which 
averages out the interaction part, leaving some room for 
controlled evolution of QC. The method is supported by 
an elegant group-theoretic framework. Here we would 
like to determine what types of reservoirs allow to ap- 
ply such method. By use of a toy model we will argue, 
that, unfortunately, typical reservoirs presented in sec. 



[II B do not allow for dynamical decoupling as proposed 
in Consider decoherence in QC with "bang-bang" 

term only, given by 



H{t) = Hcosnt 



where H is time independent, Q is large frequ 
need to derive Ya{uj). Using notation of eq. ^ 



(48) 

lency. We 
we have 

(49) 



to be weaker than the "regular" (say collisional) deco- 
herence |2j]. However, it seems that for the bang-bang 
control, the considered effect can become relevant. 

An interesting version of dynamical decoupling method 
was considered in Ref. Namely, the starting point 

was modulate the coupling constant in interaction picture 
rather than self Hamiltonian. As a result, in Schroedinger 
picture, the self-Hamiltonian was of the form uHc sinwt. 
Thus in addition to modulation, the level differences of 
the Hamiltonian were rescaled. Then the above analysis 
does not apply. However, one can show that in the cases 
where the above method works, more elementary strate- 
gies can protect the system against the reservoir. We will 
discuss it elsewhere. 



IV. EXAMPLES 

We will present here two examples. The first one will 
serve us to introduce fault map describing the action of 
noise with propagation excluded. In the second one we 
consider kicked dynamics of algorithm, and illustrate the 
difference between quantum and classical memory. 



Expanding the time dependent term into Taylor series 
and applying Fourier transform we obtain 



Ya{uj) = ^ Sa{iL>k)[coS{uj) 



ci'^{S{uj^n) + s{uj + n)) + 



(50) 



where 5 stands for (Si^'' (see eq. (p6[)), co,i are constants 
and we omitted higher harmonics. Putting Ya into for- 
mula (^^ we see that to obtain low decoherence, we need 
Raf3 (i) to be small for a; « and (ii) to have cut-off 
frequency satisfying lUc < ri. Thus one would need bell- 
shaped spectral density. The bosonic field reservoirs sat- 
isfy (i), but, as we have already mentioned, do not have 
physical cut-offs. On the other hand coUisonal or the 
coloured noise have cut-offs, yet allow for pure decoher- 
ence due to nonzero i?Q^(0). The conditions for dynam- 
ical decoupling can be in principle met in QED cavity. 
Then however, the system loses its fundamental simplic- 
ity: the quantum computer must now include modes of 
cavity, while the reservoir become the atoms building the 
cavity. 

There is another aspect of the discussed method, that 
has not been investigated so far. Namely, the time de- 
pendent Hamiltonian is obtained by use of external fields 
that can be described classically (e.g. coherent light of 
the laser beam). However, rapidly alternating field will 
get entangled with the controlled system (this is called 
quantum back-reaction) so it must be treated quantum 
mechanically. The resulting disturbance was evaluated 



A. Memory less reservoir: error map and fault map 

Memoryless reservoir has the white-noise spectral den- 
sity. Its auto-correlation function is Rap{s — t) = 
i?^^(5(s — t). Let us, for simplicity, assume independent 
interaction Raf3{uj) = RaSap- The resulting error map is 
given by 



du Sa{u, s)gSaiu, s). 



(51) 



Suppose that Sa are one-qubit operators. According to 
notation of Ref. we will call them faults - these are er- 
rors caused by the reservoir, and should be distinguished 
from the errors resulting from propagation of faults by 
algorithm (we will denote it by U{t,s) in this section). 
The propagation of the faults is described by the fact that 
the error map involves Sa in Heisenberg picture. Note, 
that the overall error is caused solely by one qubit faults 
and their propagation, averaged in time. Such evolution 
can be described by the following master equation 



— gt ^i[Ht,gt] + Lgt 



with 



L 



1 



{^*{l),g} + Mg) 



(52) 



(53) 



for ^{g) — J2a ^aQSa- Indeed, let equation (0) at be 
density matrix of QC, Ct = —iHt - where Ht is Hamilto- 
nian of algorithm and L - the above generator. Applying 
the formula (Oh in first order approximation we obtain 
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gt = U{t,to)[gt„ - -{A,gto} + 



E 



du Sa{u, s) gSa s)] 



U{t,to)[g- -Mil), q} + Mq)]- 



(54) 



Even though the memoryless case is not especially inter- 
esting from our point of view, the very form of the master 
equation (js^) is particularly useful in our context, as the 
faults themselves and the propagation are separated from 
one another. The propagation is described by the Hamil- 
tonian term, while the faults are described by the map 
$, which we will call fault map. We will derive such mas- 
ter equation in general case in next section. Here, let us 
only mention that the "fault part" of the master equa- 
tion is a sum of one-qubit operators. This means that 
environment acts locally, and the reason for occurring 
multi-qubit errors is solely due to propagation. 



B. Decoherence under kicked dynamics of QC: Error 
path interference effect 

We will assume that the gates are performed quickly, 
so that they can be generated by Hamiltonian with time 
dependence given by delta function. This will not work 
for electrodynamic or phonon vacuum reservoir, as due to 
the behaviour uj^, any rapid changes in self-Hamiltonian 
cause large or infinite contribution. Thus in this case one 
needs to work with Gaussian pulses. On the other hand, 
the kicked dynamics can be used for Lorentzian spectral 
density. 

Thus we divide computation time t into N pieces of 
length T. The Hamiltonian is given by 



Hit) 



N 

E 



Sit - ij - l)T)h, 



(55) 



where hj generates j-th step of computation. The uni- 
tary evolution is given by 



uit)^iif^,eit-3T)uij,j-i) 



(56) 



with Uij,j — 1) = e~ is the j-th step of computa- 
tion (the group of gates performed at time jr), 9 is step 
function. Consequently 



Ut 



?7(1,0) 

C/(2, !)[/(!, 0) 

UiN, N - l)UiN - 1, iV - 2) , 
...C/(1,0) 



t£ (0,t) 
t e (r, 2r) 



N 

a j,k=l 



(58) 



where S^^ = 

Rap — Ra^aP 



UjSaUj. For simplicity we assumed that 
The coefficients r° 



'jk 



are given by 



(fc-l)T 



dv Raaiu 



(59) 



For memoryless case we obtain discrete counterpart of 



N 



Ph 



(60) 



fc=i 



In both formulas we see propagation of error, this time, 
in discrete manner: 5^ denotes error resulting in propa- 
gation of fault Sa within the interval (0, jr). The differ- 
ence is that in equation ( |58| ) we have interference due to 
quantum memory. This effect of interference of error can 
remarkably modify the process of error propagation. 

The difference is similar to that between separable and 
entangled mixed states. We can view the error map as 
(subnormalized) density matrix and the Kraus operators 
as the pure components of ensemble giving rise to the 
density matrix. In general, the noise process can then 
behave as classical-like or as quantum stochastic process 
with time correlations of the form 



'^cig^ E P^^^ 



^quQ = 



E 

■ ■j»n,n, 



■i'tn)Sitj,^) 

Cin 



Sit,,)gSit,,)K 

jm I ^1 ■ ■ ■ ^m) 



■ sit.y 



Sit,J...Sit,,)gSit,,)r..Sit,y 



(61) 



respectively. In the first case we have always mixture of 
product Kraus operators, in the second one - we can have 
a map for which such representation may be impossible. 
Indeed, in equation ( |60| ) the errors occurring in different 
times do not interfere with each other, they are added as 
in classical probability calculus. 

However it seems that the crucial point is not the dif- 
ference between quantum and classical memory, but its 
range. In Sect. ^ we will exhibit remarkable connection 
between the range of the memory and occurring multi- 
qubit faults: the interference caused by long range quan- 
tum memory will give raise to nonlocal noise. 



t e {{N-1)t,Nt) 
(57) 



Substituting such dynamics into the error map of eq (25), 
we obtain 



C. Memory of electromagnetic vacuum reservoir 

Here we will calculate and discuss the memory caused 
by vacuum in free space. This is important point in our 
paper, as the vacuum is unavoidable, we cannot remove 
it. Of course we can put the QC into cavity, but then, we 
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can consider the QC plus cavity as the total system, that 
is again in free space. We then prefer the latter setup, as 
it is more fundamental. 

Consider the two-level atom and electromagnetic field 
with dipol interaction 



H = ^uJoa^ + (Ti ® dE 



(62) 



The spectral density of reservoir in temperature T is 
given by 



R 



{uj) = ^d^uj^(l+ for c^>0 

R{uj)^^d^\Lo\^( ^,.,l,_A for uj<0 



(63) 



We must now calculate autocorrelation function in zero 
temperature. Taking inverse Fourier transform we obtain 



m 



(64) 



We see that the memory scales as power of t rather than 
exponential function. Then there is no characteristic 
time. This is due to the fact, that Coulomb interac- 
tion has no characteristic range. Thus the range of vac- 
uum memory is infinite. Nevertheless, in quantum optics 
there is a notion of time of memory, which is the inverse 
of - the transition frequency for the atom, and one 
can work with Markovian master equation, within times 
much longer than I/luq- 

This does not mean that there exists some character- 
istic time scale for vacuum reservoir. Rather, in the in- 
teraction picture, the autocorrelation function acquires 
oscillating term e"^"* which effectively makes the mem- 
ory much weaker. For optical transitions, the memory 
has then almost no observable consequences. One could 
think that the memory problem is not relevant in optical 
implementations of quantum computation. However, it 
is not the case. For two atoms, there are degenerated 
levels, so that if one wants to drive gate coupling these 
levels, there is no high frequency that would diminish 
the memory effect: the frequency driving the gate and 
the frequency of the transition are the same. Then the 
time of memory and time of controlled evolution are of 
the same order. The resulting non-Markovian effects ||20| ] 
will be most relevant for many atoms (as needed in quan- 
tum computation schemes), as the number of degenerate 
levels scales exponentially with number of atoms, so that 
it will be hard to perform computation without use of 
the gates suffering from memory effects. 



V. QUANTUM LONG RANGE MEMORY 
IMPLIES NONLOCAL STRUCTURE OF NOISE 



A. Master equation 

As we have seen, to discuss the structure of noise with- 
out propagation effect, we need master equation of the 
form 



g= -i[HQc{t),g] + Ltg. 



(65) 



we will now derive such equation in second order approx- 
imation. To this end, we differentiate the equation ( [2^ ) 
obtaining 



git) = Uc{t,s)gis) + Uc{t,s) 
1 



1 



g{s)} + ^g{s) 



Uc{t, s) 



-{M,g{s)} + ^g{s) 



(66) 



We make first step of approximation by removing the 
second term. Differentiating $(t, s) and putting s = 
we obtain 



^i[HQc{t),g] -\{Tll, g} + Ttg 



(67) 



where = T(i,0) is given by 



t(i,s)(-) = ^(T„^(-)^a+^M-)Ti^) (68) 

Q/3 



with 



- Q/3 



Au Rap{u - t)Sf3iu,t) 



(69) 



On the right hand side of the equation we have den- 
sity matrices evolving according to self-Hamiltonian of 
the system. The second step of approximation is then 
to replace them with the ones subjected to full evolu- 
tion. (Better justification of such integral differential 
master equation can be obtained within cumulant ex- 
pansion method, see e.g. [^). 

Now, the decoherence is essentially characterized by 
the map Tt which we will call "fault map" in analogy 
to the error map ^t- The latter one characterizes error, 
i.e. the net effect of the "attacks" of noise (faults) and 
their propagation. The present map Tf describes solely 
the action of noise. The propagation will be due to the 
Hamiltonian term in the above master equation. 

Let us now discuss how the form of Tj depends on 
memory. The crucial term is, of course, T^^. In the case 
of no memory (and, for simplicity, independent decoher- 
ence) ,weobtain = SafsSa which is compatible with 
Sect. IV A. Assuming that each a denotes different qubit, 
we obtain that decoherence is local, and multi-qubit er- 
rors are solely due to propagation. We can then write 
the generator in the form 



(70) 
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where L" operates only on ath qubit. If we admit ex- 
ponentially decaying memory (e.g. coloured noise) with 
characteristic time comparable with the time of per- 
forming single computational step, but short in compari- 
son with time between gates, there will be some time for 
reservoir to learn something about the structure of two 
qubit Hamiltonians, so that the noise will be a sum of 
two qubit superoperators 



(71) 



Q/3 



In both cases decoherence is local. The possible faults 
are, roughly speaking, symbolized by L^,L^^ and they 
couple at most two qubits (in general, they involve so 
many qubits, as the used gates do). 

Suppose now, that the reservoir has long range mem- 
ory. Still, in each step of computation, there are two- 
qubit gates. The environment learns about their struc- 
ture a bit. The effect accumulates due to memory, so 
that after long time, the noise is strongly nonlocal: the 
generator Lt involves multiqubit operators as it keeps the 
record of all the history of self-evolution (i.e. algorithm). 
The decoherence becomes as "entangled" as the quantum 
algorithm is entangling. Note, that the effect does not de- 
pend on the time of performing gates. For short gates, 
the reservoir has short time to distinguish the temporary 
levels (driving the gate), but they are easily distinguish- 
able, as the energy differences are high. If instead the 
energy differences are low, then the time must be longer. 
This is the result of time-energy uncertainty principle. 
Thus, the relevant parameter is here "action" i.e. prod- 
uct of time of gate, and strength of applied Hamiltonian 
(energy). To perform the gate, the needed amount of 
action must be always of order of 1. Thus the general 
structure of noise does not depend very much on the way 



the gates are performed. In Sect. VI we will see that the 



strength of the decoherence does depend on it, so that 
we have the mechanism of possible optimization. 



B. Quantum memory and error correction 

Let us now discuss the implication of the above results 
on quantum error correction (QEC) method [2^j9|,p|,^ . 
The basic idea of QEC rests on two assumptions: (i) the 
environment acts locally (ii) the environment acts inde- 
pendently of self-evolution of the system. If the above as- 
sumptions are satisfied, it is indeed possible to protect the 
quantum information against decoherence by encoding it 
into highly entangled multiqubit states. In literature the 
assumption (i) was made explicitly, while the assumption 
(ii) was tacit. Thanks to assumption (i) such nonlocally 
encoded information will not be affected by decoherence 
(acting locally). The assumption (ii) ensures that the 
self-evolution will not modify the local structure of de- 
coherence. As we have argued, the latter assumption is 



never satisfied. If there is short range memory, the de- 
pendence on self-evolution is irrelevant, and QEC still 
applies. However, the interaction with vacuum, which 
is unavoidable, introduces long-range quantum memory, 
which causes environment to be rather malevolent: it 
traces self-evolution, and the more entangled the latter 
is, the more nonlocal the faults become caused by noise. 
On the other hand, the self evolution must be entangling, 
since it is to hide the quantum information just into en- 
tangled states. Thus, we have something a bit analogous 
to the Lenz rule: the structure of environment attack is 
proportional to the structure of the evolution that is to 
protect against it. It seems to be a rather pessimistic re- 
sult, as it implies, that the very idea of active protection 
against environment is fundamentally inadequate, as far 
as long quantum computation is concerned. 

However, the conclusion need not be so pessimistic. 
First of all, decoherence caused by vacuum is small in 
comparison with other sources of decoherence, e.g. de- 
phasing due to collisions (pressure decoherence), coloured 
noise or thermal noise. It seems that QEC (and hence FT 
method) can be used to fight with the latter sources of de- 
coherence, within the time window in which the vacuum- 
memory effect can be neglected. Strictly speaking other 
types of reservoirs also have long tails in memory, how- 
ever, the dominant part decays exponentially. Thus FT 
could be certainly performed within the time regime sat- 
isfying the following conditions: (i) the memory is ex- 
ponentially decaying, (ii) the memory majorize vacuum 
memory. This regime may prove sufficient to perform 
some quantum computational tasks. 



There is also another point to address |22|. The claim 
that lack of memory is sufficient for FT does not mean 
that it is in general impossible to improve FT against 
some kind of memory. As a matter of fact, the mecha- 
nism producing multi-qubit faults we described has some 
features that may be used to defend FT scheme. Namely, 
as seen from formula (|69|), the multi-qubit fault at some 
place a arise from backward propagation (with decreasing 
strength) of single qubit fault Sa ■ Now since FT method 
is able to fight with propagation forward, it cannot be 
excluded that it will be robust against faults correlated 
in such a way. However it goes beyond the scope of this 
paper. 



VI. MINIMAL DECOHERENCE MODEL FOR 
QUANTUM CONTROLLED SYSTEMS 

In this section we will postulate minimal decoherence 
model. Wc then derive formula for fidelity of quantum 
computation within the model. We shall consider a class 
of controlled open systems for which Hc{t) ^ only 
during computation of subsequent gates. The spectral 
density Rap satisfies 
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dw < oo for a some e > 0. 



(72) 



This excludes reservoirs with any pure decoherence ef- 
fects, such as gas or coloured noise. Our motivation is 
that we would like to investigate only fundamental obsta- 
cles to build quantum computer. The pure decoherence 
is not such a one, as in principle one could isolate the 
quantum computer from any influence of the above kind 
(even though it might require technology that will never 
be achieved). Since obtaining arbitrarily low tempera- 
ture is also a matter of technology, we can work with 
the only environment being vacuum. We removed time- 
independent part of self Hamiltonian, as it will only cause 
additional dissipation. 

Example of the system could be n well separated spins- 
i which interact with an electromagnetic field and pos- 
sibly with phonons. The controlled Hamiltonian Hc{t) 
can be realized by switching external magnetic fields and 
suitable interaction between spins. Another example is 
atomic degenerate metastable level, with Stark of Zee- 
man splitting used for controlling. Of course it may be 
impossible to keep possibility of controlling interactions 
(a cornerstone of a quantum computer) and remove the 
influence of pure decoherence. However, as said, we as- 
sume the most optimistic case, and keep only the obsta- 
cle that cannot be removed for fundamental reasons. By 
removing time-independent self Hamiltonian, we obtain 
that there is no decoherence on qubits that are not ac- 
tive. Thus decoherence will be caused by the vacuum 
solely due to gate operations that require to split energy 
levels during the time of performing operation. In the 
case of independent interactions of qubits with vacuum, 
this temporal difference of energy levels is made of use 
by the vacuum which forces the upper level to relax by 
spontaneous emission. In the case of collective interac- 
tion (Dicke limit) this is not the case, where one has 
decoherence-free subspaces M. 



A. Fidelity in minimal decoherence model 



We put Qc{s) = \ip){tp\ and s — 
fidelity of computation is given by 



-oo, t 



F = l-S={Uc{t, s)^\gc{t)\Uc{t, s)^j) 



-oo. The 



(73) 



We win denote Sa{t) = Sa{-oo,t) U^'^{t,~oo)Sa. To 
avoid infinite contribution from the interaction, we will 
assume that the interaction is adiabatically switched off 
at infinity (the standard method in quantum field theory 
p3[). Thus put into formula (pO|) the following operators 



Sa{t)e 



'dt. 



(74) 



Performing partial integration and then passing to the 
limit e — !■ we obtain 



oo 

lu J \at J ^ 



(co). (75) 



We then obtain the following formula for error 5 

/oo J i' 



(76) 



B. Quantum information processing in minimal 
decoherence regime 

We consider the simplest quantum gate: rotation of 
single qubit by an angle a € [0, 27r], in the case of two- 
level atom coupled to electromagnetic field in zero tem- 
perature. The time dependent Hamiltonian will be of the 
form Hc{t) = \f{t)a,^, with 



f{t)dt = a, hm fit) = 0. 

t — >±oo 



The interaction Hamiltonian is Hint = Ca: 
spectral density, as discussed earlier, is 



(77) 



and the 



R[u) 



Rquj^ for w > 
for w < 



The full unitary transformation is given by 
Uc = Uc{oo, -oo) 



e 2^- 



(78) 



(79) 



The parameter a is the amount of "action" needed to per- 
form the transformation. Roughly speaking, a is prod- 
uct of energy pumped into system during performing gate 
and the time of operation. In the case of non-rectangular 
pulse shape, the time is given by the width of pulse. Typ- 
ically, if the energy is low, then accompanied dissipation 
is small, however time of the operation is long (while 
in quantum computation, one would prefer short time 
gates). Adding more energy, we obtain short time, but 
simultaneously, enhance decoherence, causing loss of fi- 
delity. Let us now calculate the fidelity. We have 



where 



X{lu) = F^{uj)a- + F+{Lu)a^ 



/OO 
die-*"*/(t)e=^**(*) 



(80) 



(81) 



with (f>{t) = f{u)du . From the formula ( ^6[ ) one 
easily finds that any rapid variation of f{t) will result in 
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large error. Indeed, the changes produce long tails of the 
Fourier transforms -F±, which integrated with Rqu for 
Lu > give poor fidelity. In order to minimize this effect, 
we choose Gaussian shape of the pulse 



3/2 



(90) 



(82) 



where ti can be taken as the time of performing a single 
gate. Using approximate formula (j){t) ~ f{0)t + ^ one 
finds 

|i^±(u.)| =aexp{-i(..ic±a-^)'}. (83) 



Consequently the error can be estimated as follows 

S ~ i{^\a+a-\^) - Ma+\^)f)a^Rot^' (84) 

Taking the average over possible input states -ip we get 

Si - Roa^t-^ (85) 

up to a constant factor. 

We execute "algorithm" with action yl in n steps 
(gates) i.e. a = A/n ~ 0{1). It turns out that the 
formula (|7^) is additive with respect to composition of 
gates. To see it consider n gaussian pulses of width ti 
separated with time r — mti, where m » 1. This is 
smooth realization of kicked dynamics considered earlier. 
One finds 



X„(c.)=^e— ^-;7t(jr,0)ri(a 



(86) 



J-t/2 at 



(87) 



where 

r/2 

l-r/2 

Now, we have, for example 

(^■'-■'■)(t/(jr,0)V|ri^a;)ri(a;)|(7(jr,0)V'). (88) 

Taking large separation between gates, we obtain that 
the terms involving j ^ j' are averaged out by rapid 
phase rotation. The remaining term is nothing but the 
sum of single gate terms. 

Consequently, for n pulses with average error Si the 
total error is given by 



Sn ~ n6i 



(89) 



Suppose now that we would like to keep the total error 
below some threshold e then we obtain the following es- 
timation for the time of computation tc ~ Ji{m + l)ti 
(where ti is the time of single pulse) 



In conclusion, it is possible to keep high fidelity, at the 
expense of worse scaling of the physical time tc t?^"^ 
in comparison with "algorithmic" time taig ~ n. 

Finally, one should mention, that our result is not re- 
stricted to the simple one-qubit system. If one takes two 
qubits, and some two qubit gate, the reasoning will be 
similar. If the system is K-qubit, (still with, say, at most 
two-qubit gates), the only difference will be the need of 
scaling fidelity per one qubit as I/if to get high fidelity 
of the total final K-qubit state. This can be achieved 
by further slowing down the gates. Finally, to keep high 
fidelity of quantum computer in vacuum we need phys- 
ical time to scale as tc ~ n^^^^/K or tc ~ n\/v where 
V — fiK is complexity of the problem (or volume of the 
algorithm) i.e. number of steps times number of qubits 
needed to run the algorithm. 



VII. CONCLUDING REMARKS 

We have developed dynamical description of decoher- 
ing quantum computer. We have obtained dependence of 
the structure of decoherence on quantum memory of the 
reservoir. The nonlocal structure of decoherence implied 
by vacuum memory suggests, that long quantum com- 
putation at constant error rate is impossible. Also the 
high frequency "bang-bang" control is not useful for typ- 
ical reservoirs. Instead, one should deal with short time 
quantum computation by use of low frequencies. This of 
course requires optimization of the kind we performed in 
the last section. Finally, we have designed the most opti- 
mistic minimal decoherence model, and provided fidelity 
formula, in first order approximation. Using it we have 
shown how the time-energy trade-off influences scaling of 
physical time of computation. 

In general we have argued that in description of quan- 
tum computing the Markovian approach fails, as far as 
linear coupling to boson field is considered. The ef- 
fects are relatively small for coupling to electromagnetic 
vacuum (though certainly much stronger than the non- 
observable non-exponential decay of radiating atom). 
Our results can have practical meaning for the systems 
interacting with phonons, like e.g. quantum dots (cf. 
p6| ), where the coupling is much stronger. Then the 
memory effects can be of similar order as other sources 
of decoherence. There are also less fundamental but prac- 
tical sources of memory such as e.g. heating mechanism 
in ion traps [^ . 

Finally, we believe that our non-Markovian dynam- 
ical description of quantum computing will have some 
meaning for future practically useful implementations of 
quantum computer involving large number of qubits and 
relatively long time of computation. 
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